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Abstract. Graphite is a layered material showing a strong anisotropy. Among the 
unconventional properties reported by experiments, the electronic transport along 
the c-axis, which has direct implications in order to build graphitic devices, re- 
mains a controversial topic. We study the influence of inelastic scattering on the 
electron tunnelling between layers. In the presence of electron electron interac- 
tions, tunnelling processes are modified by inelastic scattering events. 



1 Introduction 

Graphite is composed of stacked layers of two-dimensional graphene. It is known that layered 
materials have been object of intensive study since they present important physics. A key 
feature of many layered materials is the anisotropy exhibited by its transport properties: while 
being metallic within the layers, the transport in the c-axis, perpendicular to the layers, may be 
coherent or incoherent and undergo a crossover with temperature from one regime to the other, 
thus changing the effective dimensionality of the system |1|2) . Even when coherent electron 
excitations can be assumed within individual layers, there is no consensus about over what 
length and time scales the excitations are coherent between layers pQ. Electron correlations 
play as well a key role in the physics of layered materials since correlation effects increase as 
dimensionality decreases, therefore dimensionality is crucial for the electronic properties and 
to choose the appropriate model to study the system. Unusual properties are derived from 
the anisotropy and periodicity along the axis perpendicular to the planes i.e. the structure of 
collective excitations absent in two dimensional (2D) and three dimensional (3D) electron gases 
[3] These unconventional properties are distinctly different from that of the traditional Fermi 
liquids. 

Graphite, an hexagonal layered material, presents an intraplane hopping much larger than 
the interplane hybridization. Many of the transport properties established in the past for this 
well known material are being questioned at present. Recent conductivity measurements reveal 
a suppression of the c-axis conductivity much larger than what would be predicted by the band 
calculations of the interlayer hopping [4] . Band structure calculations are also challenged by very 
recent claims of observation of quantum Hall plateaus in pure graphite [5] . The unconventional 
transport properties of graphite such as the linear increase with energy of the inverse lifetime 
[6|7j . suggest deviations from the conventional Fermi liquid behavior, which could be due to 
strong Coulomb interactions unscreened because of the lack of states at the Fermi level [8 9 . The 
experimental isolation of single graphite layers (graphene) [lOj has enhanced the importance 
of graphite opening new possibilities in the application field. The properties of few of single 
layers of graphite have been measured and the experiments confirm the theoretical predictions 
of a physics governed by the two-dimensional Dirac equation and a linear dispersion in the 
proximity of the K points [1011 lj 
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By assuming that electron correlations modify the in-plane electron propagators |12j . we 
show that even in the clean limit, many-body effects can alter the coherent contribution to the 
out of plane electron hopping. The clean limit is defined as that in which the length scale over 
which electrons remain coherent within the layers, diverges. It is shown that, for certain models 
of correlated electrons [819113] , the interplane hopping between extended states can be a relevant 
or an irrelevant variable, in the renormalization group (RG) sense, depending on the strength of 
the coupling constant. The interlayer coupling in graphite is still under debate. Band structure 
calculations predicted an interlayer coupling much larger than the values deduced by c-axis 
conductivity measurements [14115] . Furthermore, the traditional assumption of graphite as a 
semimetal does not consider the electron-electron interactions which, due to the low density 
of states, remains unscreened [14I16J . The scheme used here is based on the RG analysis as 
applied to models of interacting electrons [17] . 

Two-dimensional systems that support low-energy excitations which can be described by 
Dirac fermions have recently attired a renewed interest. Besides quasiparticles of planar zero- 
gap semiconductors as graphite [8|9] there are examples as the flux-phase of planar magnets 
[18j , nodal quasiparticles in d-wave superconductors [19)20] , or the insulating spin density wave 
phase of high T c superconductors [21] . 

Long range Coulomb interactions and disorder are modelled in these systems by means of 
gauge fields coupled to the Dirac fermions [22] . Most of these systems show anomalous transport 
properties ranging from mild departures of Fermi liquid behavior as in pure graphite, to the 
total destruction of the quasiparticle pole. Disorder can be modelled as random gauge fields 
coupled to the Dirac quasiparticles. It affects the quasiparticle Green's function in a computable 
way what in turn modifies the interlayer tunnelling. 



2 The method of calculation 

We analyze here the interlayer hopping in a layered material such as graphite. In the pres- 
ence of electron-electron interactions, tunnelling processes are modified by inelastic scatter- 
ing events. The influence of inelastic scattering on electron tunnelling has been studied, us- 
ing equivalent methods, in mesoscopic devices which show Coulomb blockade [23], Luttinger 
liquids 24 25 26 27J, and dirty metals[28 . The simplest formulation of the method replaces the 
excitations of the system (such as electron-hole pairs) by a bath of harmonic oscillators with 
the same excitation spectrum. This approach can be justified rigorously in one dimension, and 
is always an accurate description of the response of the system when the coupling of the quasi- 
particles to each individual excitation is weak although the net effect of the environment on 
the system under study can be large The expression for the coupling between the electrons 
and the oscillators is obtained by assuming that the oscillators describe the charge oscillations 
of the system. Then the squared coupling is related, using perturbation theory, to the charge- 
charge correlations in the electron gas. This is consistent with the assumption that the modes 
in the environment are weakly perturbed by their interaction to the low energy electrons whose 
tunnelling properties are been considered: 



where c](ci) creates (destroys) an electron at site i, and /5k describes the charge fluctuations 
of the environment, which are to be described as a set of harmonic modes as stated above, 
he interaction in cq.JT]) is the simplest, and most common, coupling between the tunnelling 
electrons and the excitations of the system, which is spin independent. Other spin-dependent 
or more complicated couplings can also be taken into account, provided that the appropriate 
response function is used. 

Since the excitations of the system (electron-hole pairs, plasmons) are modelled as bosonic 
modes, one can write an effective electron-boson hamiltonian of the type: 




(1) 



k 
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= ^2 t io4 c i +^u k b\b k + ^2g k , l 4 c i( b l + b k) (2) 

where H e i ec describes the individual quasiparticles, H env stands for the set of harmonic os- 
cillators which describe the environment, and Hint defines the (linear) coupling between the 
two. 

The Ouipk) are boson creation(destruction) operators, the t,j describe the electronic hop- 
ping processes. The information about the interaction between the electron in state i and the 
environment is defined by the spectral function [29] 

Ji(u>) = 2J l5fc,i| 2 <5(w - cj fe ). 
k 

Using second order perturbation theory and eq.([T]), we can write |23t29] : 

JiH = ^^ 2 (k)Im X (k, W ) (3) 

k 

where x(k, to) is the Fourier transform of the density-density response of the system, (/Sk(i)/>-k(0)). 
The electron-boson interaction leads to a Franck-Condon factor which reduces the effective 
tunnelling rate. The Franck-Condon factor depends exponentially on the coupling between the 
particles and the oscillators, and when it diverges additional self consistency requirements be- 
tween the hopping amplitudes and the oscillator frequencies included in the calculation have 
to be imposed. The electron propagators acquire an anomalous time, or energy, dependence, 
that can be calculated to all orders if the state i is localized, or, which is equivalent, neglecting 
the hopping terms in eq. ([2|) ■ The present scheme can be considered a generalization to layered 
systems of the approach presented in reference[25] for a Luttinger liquid, where the scaling of 
tunnelling is expressed in terms of an effective density of states. To second order in perturbation 
theory, the electron propagators in a localized system are 



(4(t)Ci{t')) ~ (4(tMt'))o x 
exp < — / duj 



e iu(t-t') 



(4) 



where {c\{t) Cl {t')) Q ~ e l£ * (*-*') is the Green's function in the absence of the interaction. In 
contrast to previous work, we analyze tunnelling between coherent extended states. In order to 
do so, we need to generalize eq.Q to this case. We then need to know the Green's function of 
coherent states in the individual layers, G(k, w), including the correction due to the interaction 
to the environment. We firstly assume that eq.((4]) also holds in a system with extended states. 
For a standard metallic system, we must insert (cf (t)cj(£'))o ~ l/(t — f') in eq.Q. It can be 
shown that this approximation is exact at short times, W <C (t — t')^ 1 <C A, where W is an 
energy scale related to the dynamics of the electrons, and A is the upper cutoff in the spectrum 
of the environment. This expression can be generalized, taking into account the spatial structure 
of the coupling to: 



t-t> 



xexpj / / d 2 rdVd 2 ke lk ( r - r '> 



l-e % 



j(t-t') 



(5) 



where J? is the region of overlap of the wavefunctions !Pj(r) and #j(r). This expression, which can 
be seen as the exponential of the leading frequency dependent self-energy correction to the elec- 
tron propagator, has been extensively used in studies of tunnelling in zero dimensional systems 
(single electron transistors) which show Coulomb blockade [23]. one dimensional conductors [27]. 
and disordered systems in arbitrary dimensions [28] . 
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The effective interaction can, in turn, be written in terms of the response function as: 

K//(k,^) = ^ 2 (k)Im X (k, W ) (6) 

The time dependence in eq. ([5]) is determined by the low energy limit of the response 
function. In a Fermi liquid, we have 

Imx(k,w) w a(k)|w| u < E F , (7) 

where Ep is the Fermi energy. We then get: 

lim (&WV))~ T* !m +a ( g ) 
(t-t')— >oo (r — V ) i+Q 

where 

a = f f f d 2 rd 2 r'd 2 ke lk{r - r ' ) V 2 {k)a(k) (9) 

the parameter a gives the correction to the scaling properties of the Green's functions. Inte- 
gration in k is restricted to |k| <C L^ 1 where L is the scale of the region where the tunnelling 
process takes place. The value of L is limited by the length over which the phase of the electronic 
wavefunctions within the layers is well defined. 

We can now use eq.® to analyze the interlayer tunnelling by applying renormalization 
group methods. The simplest case where this procedure has been used is for the problem of an 
electron tunnelling between two states, i and j, which has been intensively studied 30 31J. We 
integrate out the high energy bosons, with energies A—dA < ujk < A and rescaled hopping terms 
are defined. As mentioned earlier, eq.([8]) is valid for this range of energies. The renormalization 
of the hoppings is such that the properties of the effective Hamiltonian at energies ui <C A 
remain invariant. If the hoppings ty are small, any physical quantity which depends on them 
can be expanded, using time dependent perturbation theory, in powers of: 

t^fteWcJtOcitO) « t?.(ct(f)c-(t')>< Ci (t)c)(t')) (10) 



The integration of the high energy modes implies that the terms in eq. (|10|) are restricted to 
t < A^ 1 , or, alternatively, the time unit have to be rescaled[32], r' = re dA ^ A , where r ~ A^ 1 . 
Using eq.®, the condition of keeping the perturbation expansion in powers of the terms in 
eq. (|10[) invariant implies that: 

t%^t%e^ 2 ^ (11) 
which can also be used to define the scaling dimension of the hopping terms. Finally, 

aft^M) = _ a t_^ (12) 

where I = \og(Ao/A), and Aq is the initial value of the cutoff. 

This approach has been successfully used to describe inelastic tunnelling in different situa- 
tions in [12|23|24|25|26|27|28j . 

The analysis which leads to eq. (fT2")) can be generalized to study hopping between extended 
states, provided that we can estimate the long time behavior of the Green's function, 

G(k,t-t') = (ct(t)c k (t')). (13) 

We assume that, in a translationally invariant system, there is no dependence on the position 
of the local orbital, i. This result implies that the frequency dependence of the Green's function, 
in a continuum description, can be written as: 

lim G(r-r',w) oc U Q (14) 

|r-r'|->0 
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Equation (|T3|) is related to eg. (fl4| . by: 

lim G(r -r',uj) = [ d D kG(k,Lu) (15) 

|r-r'|->0 J 

where D is the spatial dimension. In the cases discussed below, the interaction is instantaneous 
in time, and the non interacting Green's function can be written as: 

G (k,u) oc -T ( (16) 
uj \ u> J 

where z = 1,2 depending on the dispersion relation of the given model. In the following, we 
assume that the interacting Green's function has the same scaling properties, with the factor 
uj^ 1 replaced by lu~ s in eq. (|16p . where 5 depends on the interactions. This can be shown to be 
correct in perturbation theory to all orders, in the models studied below, which describe the 
physics in the proximity of critical points, because the corrections depend logarithmically on 
u> (it is a well known fact for the Luttinger liquid). Then, using eqs. (fT4|) . (fT5|) and (fT6|) . we 
obtain: 

G(k,uj) oc \w\ a - D '*r(£) (17) 



and T(u) is finite. Thus, from the knowledge of the real space Green's function, using eq.(|l]), 
we obtain a, which, in turn, determines the exponent a + D/z which characterizes G(k,u>). 
Generically, we can write: 

G loc , ext {u>) ~ M<W* t (18) 

where the subindices ioc, ext stand for localized and extended wavefunctions. In terms of these 
exponents, we can generalize cq. (fT"2"|) to tunnelling between general states to: 

— —oi oc ^ ext — l^iyj 

Before proceeding to calculations of Si oc and S ex t for various models, it is interesting to note 
that, in general, the response function of an electron gas in dimension D > 1 behaves as 

lim x( k >^) ~ M/|k|, 
w— >0, k — >Q 



so that, from eq.< 



lim a ~ L^-V 

L — >oo 



Thus, for D > 1, the contribution of the inelastic processes to the renormalization of the 
tunnelling vanishes for delocalized states, L — ► oo. 

It is easy to show that, in an isotropic Fermi liquid in D dimensions, linii^oo a oc L X ~ D , 
where L is the linear dimension of the (localized) electronic wavefunctions tf'(r). This result is 
due to the dependence on k of the response functions which goes as 

Im%(k, ui) 



k D F ^\k\ 



Thus, for D > 1, we recover coherent tunnelling in the limit of delocalized wavefunctions. 
In one dimension, one can use the non interacting expression for hnxo(k,u)), to obtain: 



d 2 rd 2 r'd 2 ke lk ^- r "> dull- e 1 ^'^ 



LU 2 

I U V JO t-t' «i/«F , on ^ 

K \e;) X \ log[v F (t - t')/L] t - t> > L/vf (M) 
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where we have assumed a smooth short range interaction, parametrized by U . Hence, the 
Green's functions have a non trivial power dependence on time, even in the L — > oo limit, 
in agreement with well known results for Luttingcr liquids 27 . In order to obtain the energy 
dependence of the effective tunnelling between k states near the Fermi surface, one needs to 
perform an additional integration over dr. In general, near a scale invariant fixed point, uj oc |k| z , 
and for a ID conductor one knows that z — 1. Hence, 

Im G(uj, k F ) oc uj- z+a ~ uj- 1+a . 

The flow of the hopping terms under a renormalization group scaling of the cutoff is [26] : 

d(t/A) _ ( — a localized hopping 

dl 11 — a extended hopping ' 

where t denotes a hopping term, between localized or extended states. In the latter case, the 
hopping becomes an irrelevant variable [24|33|34] for a > 1. 



3 Graphite layers 

The simplest two dimensional model for interacting electrons where it can be rigourously shown 
that the couplings acquire logarithmic corrections in perturbation theory is a system of Dirac 
fermions (e^ = v,p|k|), with Coulomb, l/\r — r'|, interaction. This model can be used to describe 
isolated graphene planes [8,9], formed by carbon atoms arranged i the honeycomb lattice, and 
can help to understand the anomalous behavior of graphite observed in some experiments |6l7j . 

In order to apply the procedure outlined in the previous section, one needs the Fourier 
transform of the interaction, 

K//(k) = e 7(eo|k|), 

where e is the electronic charge, and eo is the dielectric constant, and the susceptibility of the 
electron gas. For a single graphene plane, this quantity has been computed in [Sj and is: 

Xo (k, t j) = l- r . |k|2 (22) 

These expressions need to be inserted in equations <j6j> and ([5]). Alternatively, we can use the 
RPA, and include the effects of interplane screening: 



XRPA{k,u) = 

sinh(|k|d) 



[cosh(|k|d)+f^siah(|k|d) X o(k,w)] 



(23) 



where d is the interplane spacing. The imaginary part, Im%o(k, uj), is different from zero if 
uj > uj?|k|. 

For simplicity, we consider the expression in eq. (!22j) . as it allows us to obtain analytical 
results. We cut off the spatial integrals at a scale, L, of the order of the electronic wavefunctions 
involved in the tunnelling. Performing the same computation as in the case of D=l we obtain 
an expression similar to that in eq. ([20|) except that the prefactor (U/Ep) 2 is replaced by the 
squared of effective coupling constant of the model, e 4 /(egvp) 2 . Thus, also in the graphene 
model, the propagators acquire an anomalous dimension what was advocated in [8] as pointing 
out to a departure of the model from Fermi liquid behavior. As in ID, the value of the exponent 
z which relates length and time scales is z = 1. The scaling of the hoppings now are: 

A d{t/A) _ ( -1 - a localized hopping . . 

dA 1 1-Q extended hopping ^ ' 



The extra constant in the first equation with respect to eq. (j2Tj) reflects the vanishing of the 
density of states at the Fermi level for two dimensional electrons with a Dirac dispersion relation. 
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In graphite, the dimensionless coupling constant, e 2 /vF, is of order unity. Under renor- 
malization, it flows towards zero[8]. Thus, despite the departure from Fermi liquid behavior, 
interplane tunnelling is a relevant variable and a coherent out-of-plane transport in pure clean 
graphene samples should be observed 2 2J. 



3.1 Influence of disorder. 

In previous section, the calculation has been done in the clean limit. This picture can change 
in the presence of disorder that is known to change the anomalous dimension of the fields. 
Disorder can be incorporated in the Fermi liquid two dimensional systems in the clean limit. 
In the renormalization group scheme disorder can be included by the introduction of random 
gauge fields 35 36 1. In fact, this is a standard procedure in the study of the states described by 
the two-dimensional Dirac equation associated to random lattices or to integer quantum Hall 
transitions and, usually, the density of states at low energies is increased. Topological disorder, 
as pentagons and heptagons in the hexagonal lattice, have been observed in graphite and 
carbon nanotubes samples. A random distribution of topological defects can be described by a 
non-Abelian random gauge in the calculation[36j and this features will influence the interlayer 
transport properties |34|. Disorder in general, will affect the electronic properties in different 
ways. Ferromagnetic behavior of graphite samples, enhanced by proton irradiation, has been 
reported [37]. The formation of vacancies, cracks or voids in the graphite lattice, due to the 
proton bombardment, favors the magnetic order in the sample due to the localized states 
which arise at the borders of the defects [35] , where the change of coordination of the carbon 
atoms gives rise to the existence of unpaired spins. The incompletely screened electron-electron 
interactions will induce the ferromagnetic correlation between these unpaired spins [38 39 in the 
layers. The influence of the long-range Coulomb correlation in the interlayer coupling remains 
to be completely understood. Recently The possibility of making atomic scale devices based on 
graphene layers has been suggested by the experiment. The control of the electronic structure 
of a bilayer graphene has been reported [40]. The gap between the valence and conduction bands 
has been controlled by selectively adjusting the carrier concentration of each layer, to change the 
Coulomb potential, at higher electron density the overlap between orbitals of adjacent layers 
increases, thus increasing the interlayer coupling. The interlayer interaction varies then as a 
function of the electron concentration in the layers. 

The possibility of controlling as well the spin transport in graphene and in CNTs, is being 
currently investigated. From the theoretical point of view, the inclusion of the spin-orbit inter- 
action in the Hamiltonian|4I 42 43 may explain some results of long spin correlation length in 
CNTs contacted to ferromagnetic leads. . 

4 Conclusions 

The potential applications of graphite sheets together with those of CNTs have caused a re- 
newed interest in carbon-based materials. In order to understand their behavior under external 
probes, it is important to know the dielectric response, which is related to the collective ex- 
citations of the electrons which combine in-plane and inter-plane interactions. Therefore the 
exact nature of the interlayer coupling is needed. We have analyzed the relation between the 
interlayer tunnelling and the inelastic processes in the layers, in the clean limit at low en- 
ergy. Our results suggest that, when perturbation theory for the in-plane interactions leads to 
logarithmic divergences, the out-of-plane tunnelling acquires a non trivial energy dependence. 
This anomalous scaling of the interlayer hopping can make it irrelevant, at low energies, if 
the in-plane interactions are sufficiently strong[34]. A relation between the nature of in-plane 
excitations and the perpendicular transport has been reported from ARPES and transport 
measurements in strongly correlated layered materials: when coherent quasiparticles form in 
the planes, the perpendicular transport becomes metallicpQ. In graphite, we find that, despite 
the deviation from Fermi liquid behavior, coherent out-of-plane transport can occur. This re- 
sult makes contact with experiments which suggest that the temperature dependence of the 
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out-of-plane resistivity is directly correlated to that of in-plane resistivity in graphite [44]. As 
well, the monotonic increase of the interlayer interaction as a function of electron concentration 
in both layers recently observed in a bilayer graphene|40j agrees with our results. 
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